Exciton resonances quench the photoluminescence of zigzag carbon nanotubes 
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We show that the photoluminescence intensity of single-walled carbon nanotubes is much stronger in tubes 
with large chiral angles - armchair tubes - because exciton resonances make the luminescence of zigzag tubes 
intrinsically weak. This exciton-exciton resonance depends on the electronic structure of the tubes and is found 
more often in nanotubes of the + 1 family. Armchair tubes do not necessarily grow preferentially with present 
growth techniques; they just have stronger luminescence. Our analysis allows to normalize photoluminescence 
intensities and find the abundance of nanotube chiralities in macroscopic samples. 
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A major challenge in research on single walled carbon nan- 
otubes (SWNT) is to control and measure the nanotube chiral 
indices on bulk, macroscopic samples. The chiral index (n,m) 
fixes the nanotube's diameter and chiral angle. These two pa- 
rameters determine all properties of a nanotube, in particu- 
lar its electronic structure. 1 1 ] Photoluminescence (PL) from 
single-walled carbon nanotubes in solution decreases strongly 
in intensity for nanotubes with small chiral angles (zigzag 
tubes) 0]. This has been interpreted as reflecting the 
abundance of (n,m) nanotubes[2]. If correct, the interpreta- 
tion has far reaching consequences, because it implies that 
all present growth techniques strongly favor armchair over 
zigzag tubes. Is, however, the luminescence cross section re- 
ally independent of the chiral angle? Can we expect constant 
maximum luminescence intensities when comparing two nan- 
otubes of different chirality? 

In this paper we show that luminescence strongly favors 
a subset of nanotubes — tubes with a small ratio between 
their second and first transition energy. This, in particu- 
lar, implies a stronger luminescence intensity for large chi- 
ral angles and nanotubes from the —1 family as observed 
experimentally. [2, 4] The chirality dependence arises from an 
exciton-exciton resonance. It also shifts the experimental op- 
tical transition energies to the red when compared to the bare 
exciton energies. From the experimental data we calculate 
the exciton energies and maximum luminescence intensities 
of more than 40 nanotube types. 

We consider first luminescence in two specific nanotubes 
(11,1) and (10,2) with similar diameter but very different ex- 
citon behavior. Our argument is based on a key observation. 
The optical transition energies are often represented by the 
'Kataura plot' in which the energies of each subband vary 
roughly inversely with nanotube diameter. The transition en- 
ergies deviate systematically above and below this trend, ac- 
cording to families of chiral indices, ill We show that this 
deviation leads to an extra exciton decay channel in tubes with 
small band gaps. 

Figure E a ) shows the PL process in SWNTs. A photon 
hV22 creates an exciton eh22 in the second subband of the tube, 
where the index 2 refers to subband 2. The exciton relaxes to 
the lower subband state ehn and recombines to the ground 
state, emitting the photon hv\\. Strong PL occurs if I1V22 cor- 



responds to a singularity in the excitonic density of states, see 
Fig. Ed). 

FigureE a ) changes fundamentally if we allow the presence 
of two excitons in subband 1, which we denote by 2ehn.|@] 
In a (10,2) nanotube, the eh22 energy £22 is unusually low 
so the 2-exciton state lies above the eh22 state. The standard 
picture is retained, Fig.E D X an d the exciton just decays into 
subband 1. On the other hand, in the (11,1) nanotube, eh22 
is high in energy and 2ehn lies below it. The eh22 exciton 
decays into the ehn state with energy E\\ and thereby liber- 
ates enough energy to create a second exciton ehn, Fig. [He). 
There are two crucial points about the eh22 — * 2ehn decay 
channel: Whether it is allowed energetically or not depends 
on the nanotube structure. Two seemingly similar nanotubes 
can show very different exciton dynamics. Second, the 2eh22 
state has a singular energy dependence [Fig. Ed)]; therefore 
the higher-order process strongly affects the nanotube optical 
properties. 

In bulk SWNT samples the PL of (1 1,1) nanotubes is found 
to be three times weaker than from (10,2) nanotubes, see 
Fig. |2ja). This factor was observed in nanotubes grown by 
different methods|2, 4j], so it is unlikely to be due to a chiral- 
ity or diameter dependence in the growth process. The (11,1) 
and (10,2) nanotubes have very similar diameters (9.0 and 8.7 
A) and similar, small chiral angles (4.3° and 8.9°). They are 
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FIG. 1: (Color online) (a) Standard picture of eh22 excitation from 
the ground state (GS) by photon absorption (full arrow), exciton re- 
laxation to ehn (broken), and emission (full) of h\\\. Carbon nan- 
otubes fall into those (b) with only one decay channel and those (c) 
where eh22 can decay into two ehn; (d) excitonic density of states 
(DOS). 
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FIG. 2: (Color online) (a) Columns: excitation wavelength with the 
strongest PL and the intensity of tubes in bulk samples. 1 2, 4] Lines: 
calculated profiles (see text), (b) Green, italic: p = 0, red, bold: +1, 
and blue: —1. 



neighbors on the standard (n,m) plot of a graphene sheet in 
Fig. |2t)). However, one characteristic clearly distinguishes 
(10,2) and (11,1) tubes - they belong to different index fami- 
lies. 

SWNTs are characterized by a 'family index' p = (n — 
m)mod3.|5] Tubes with p = are metallic. |1] Tubes with 
p = +1 or —1 are semiconductors. For +1 tubes like (11,1), 
the eh22 transition energies £22 lie above the averaged Kataura 
trend, for —1 tubes they lie below this trend, Fig. Ufa). 
This arises from the trigonal warping of the graphene band 
structure. |5| |8| This difference carries over into the exciton 
energies. Tubes with p = +1, like (11,1), have £22 > 2£n 
as in Fig. [He), while those with p = —1, like (10,2), have 
£22 < 2£n as in Fig. [0b). Now, nanotubes in the +1 fam- 
ily have the extra eh22 — > 2ehn relaxation channel, which is 
forbidden in —1 tubes. We argue that this weakens the maxi- 
mum luminescence intensity in +1 tubes by broadening their 
absorption line-width. 

To test this idea, we calculated the absorption spectrum for 
the efi22 exciton by a Green function method. 0] 
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with a dielectric screening a — 0.150 H, 
coupling A a 2 E\\ 0.01 eV 2 , and an eh22 — 
B 2 a 4 = 0.1 eV 2 



an eh 2 2 —►ehn 
2ehn coupling 
The parameters were adjusted to fit the 
electron-hole decay times in single-walled carbon nanotubes 
and graphite. 11(1 HHl En (£22) is the energy of the ehn (eh22) 
exciton. 

Figure|3a) shows the calculated absorption as a function of 
excitation wavelength. The (10,2) absorption shows a single, 
slightly red-shifted Lorentzian (see arrow at £22)- However, 
in the (11,1) nanotube, the decay to 2ehi 1 changes the absorp- 
tion; the strongly red-shifted 640 nm peak has a weak side- 
band at higher energies. Identifying the two narrow peaks in 
the calculated spectrum of the ( 1 0,2) and (11,1) tube with their 
measured PL excitation energy and intensity, we get excellent 
agreement between theory and experiment, Fig. |2ja). We find 
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FIG. 3: (Color online) (a) Closed symbols show the bare exciton 
energies; open symbols were calculated. Full lines connect tubes 
in a branch (2n + m =const.), dashed lines tubes with largest chiral 
angles (close to armchair); squares are zigzag tubes, (b) Ratio of £22 
and E\ 1 for three branches. The dashed blue line is the as measured 
ratio IiV22/hVu. 



that the maximum PL intensity of a (1 1,1) tube is intrinsically 
weaker than that of a (10,2) nanotube. 1 121 

We now generalize our findings to arbitrary carbon nan- 
otubes. First we calculate the bare eh22 energies rather than 
the renormalized energies I1V22 measured by optical spec- 
troscopy. For £n we use the measured transition energies, 
since this exciton is not renormalized. ||2j |3[ £22 we find in 
a self-consistent routine by requiring the maximum absorp- 
tion probability in Eq. Q to occur at hV22- I1V22 is, in general, 
smaller than £22 because of the renormalization, see Fig.Qand 
the supplementary material. Additionally, we calculated 
£22 within the third-order tight-binding approximation;! 14] 
we used the parametrization by Kane and Mele to account for 
electron-electron and electron-hole interactional. 

Figure |3 a) shows the bare exciton energies. Above the 
isotropic line (black dashed), the agreement between theory 
and experiment is excellent; below there are deviations of 10 
- 20 % that arise from curvature. The agreement between the- 
ory and experiment systematically improves for large diam- 
eters and chiral angles, see Fig. |3b). For close-to-armchair 
tubes (dashed lines) experiment matches theory above 1 1 A. 
In contrast, the as-measured hV2i/hV\\ ratio is constant and 
10 % smaller than the exciton ratio, which is called the 'ratio 
problem'. J2H 

The PL intensity depends on the product of the absorption, 
relaxation, and emission probability. For a specific tube the 
PL intensity as a function of excitation energy simply follows 
Eq. Q. For different tubes we assume constant thermaliza- 
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FIG. 4: (Color online) (a) Calculated photoluminescence as a func- 
tion of excitation and emission wavelength. The dashed white 
line shows the near-armchair direction, (b) Measured (Fig. 3(d) in 
Ref. 1 4]) and calcuated PL intensity for HiPCo samples as a function 
of nanotube diameter and chiral angle. 



tion rates. The optical matrix elements depend on diameter as 
1 /d; their dependences on chiral angle and family cancel, be- 
cause of opposite trends for the second (absorption step) and 
first (emission) subband. The luminescence intensity hence 
follows the absorption probability [Eq. Q] weighted by l/d 2 
and the abundance of nanotube chiralities. Here and in the 
remainder of the paper we consider a simulated sample of 78 
nanotubes types with mean diameter of 9.5A, a Gaussian di- 
ameter distribution of width of 2A, typical of HiPCo samples, 
and no preferred chirality. 

The calculated luminescence map in Fig. @Ia) agrees well 
with the experimental false color plot of Bachilo et al.\2]. 
For comparison we provide a PL map for the same ensemble 
where we neglected the eli22 — > 2ehi i decay [B = in Eq. Q] 
as a supplement. 1 13] In Fig. |4ja) the intensity is strong for 
large chiral angles in both semiconducting families. These are 
the PL peaks close to the armchair direction, see dashed line in 
Fig- Eta). Above this line is the —1 family of semiconducting 
tubes. There are 15 peaks clearly visible in this region. Below 
the dashed line the plot looks very different in both theory and 
experiment from above the line. Although there are 14 tubes 
in this region, we predict only four well defined peaks (com- 
pare supplementl 13]). They are +1 tubes with large chiral 
angles. For smaller chiral angles, trigonal warping lowers the 
ehn energy, but raises eh22- The eli22 — > 2ehn resonance sets 



> 


1 .4 


o 




LU 


1.2 


c= 


o 




£ 


1.0 


'o 




X 




LU 






0.8 



1 100 

10 

gov 

.0.05 
■ <0.01 



200 250 300 350 400 

RBM frequency (cm 1 ) 

FIG. 5: (Color online) Calculated intensity of the RBM in reso- 
nance with ehn for HiPCo nanotubes. The white line shows the 
near-armchair direction. 



in; it blurs and broadens the absorption spectra, see Fig.^a). 
In the luminescence map this creates the vertical streaks. 

Figure|4lb) compares the measured and calculated lumines- 
cence intensities; we see excellent agreement. The intensity 
of —1 tubes (blue) is almost independent of chiral angle. In 
contrast, +1 tubes (red) show strong luminescence for large 
angle, but are very weak towards the zigzag direction. In ad- 
dition, luminescence hardly sees some semiconducting nan- 
otubes at all. For example, the luminescence of a (13,0) tube 
(small chiral angle) is 5 times weaker than for a (7,5) tube. 
Thus, luminescence is strongly biased towards large chiral an- 
gles and nanotubes with (n — m) mod 3 = — 1 . 

For one nanotube, the (11,0) tube, we predicted a strong 
luminescence intensity although it is absent in the measured 
spectra [Fig.0Jb)]. This could imply a small number of (1 1,0) 
nanotubes in the sample. Another explanation, however, are 
optically inactive excitons below ehi 1 . fl5Tl . The (1 1,0) tube is 
a singular case in the experimental data as well. The two other 
— 1 nanotubes with very small chiral angles [( 1 0,2) and (12,1)] 
have strong intensities. Two other points are noteworthy: first, 
the calculations do not predicted a constant background above 
1000 nm emission wavelength. Further studies are desirable 
to clarify the experimental background. Second, in our cal- 
culation, there are no features below 400 and above 850 nm 
excitation wavelength, because we considered only the eh22 
and the ehn state and their interactions. 

How can we further verify the model experimentally and 
what are the practical implications for finding nanotube abun- 
dances from optical spectroscopy? The most rigorous test is 
to compare PL intensities with a chirality distribution obtained 
from a non-optical technique such as electron diffraction. This 
would establish an experimental normalization for the lumi- 
nescence intensities in addition to the theoretical factors given 
by us. lll 311 Time-resolved spectroscopy could be used to ob- 
serve the two distinct eli22 — > ehn an d en 22 — * 2ehn decay 
channels. The challenge is the weak luminescence signal in 
tubes with the latter decay process. Another prediction from 
our model is a strong difference between the Raman cross sec- 
tion in resonance with eli22 and ehn. 

The Raman cross section is proportional to the square of 
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the absorption strength. Raman scattering in resonance with 
eh22 therefore shows a similar depedence on tube type as pho- 
toluminescence. It is not exactly the same, because of the 
squared absorption probability and electron-phonon coupling, 
see Ref. 11611 . For Raman scattering in resonance with ehn, 
however, we predict a remarkably straightforward way to ex- 
tract the chirality abundance, because there are no excitonic 
states below the first subband exciton. 

We calculated the Raman spectra of the radial breathing 
mode (RBM) in resonance with ehn from non-orthogonal 
tight-binding. B 1 711 Although this calculation neglects electron- 
hole interaction, the dependence of the optical and electron- 
phonon matrix elements on diameter and chiral angle should 
be well described by a one-electron model. The Raman in- 
tensity map is shown in Fig. [5] In contrast to luminescence 
[Fig.Ua)], Raman scattering shows well resolved peaks for 
both semiconducting families and all chiral angles. 

The maximum Raman intensity in Fig.[5]varies for different 
RBM peaks. This comes from the Gaussian diameter distri- 
bution and the dependence of the Raman matrix elements on 
chiralitv. 1 1711 The latter is described excellently by simple an- 
alytic functions of the nanotube diameter d, chiral angle ©, 
and family p. The square root of the Raman intensity follows 
to very good approximation (less than 10 % deviation) 



X= l+p4.82-10 _ 



4.59- \0- 2 d 



nm 



1 + (2.66 -10 
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+ /J6.18- 10^)0+ 1.06- 10^© z Xo, (2) 

where is given in degrees and Xo is constant. Infrared Ra- 
man scattering could thus verify our model of exciton decay. 
Additionally, the ehn -resonance intensities could be used to 
normalize PL intensities and their dependence on family and 
chiral angle. 

In conclusion, luminescence has a systematically higher 
cross-section for SWNTs with large chiral angles and tubes 
from the —1 family. This arises from an additional decay 
channel when the exciton of the second subband has more 
than twice the energy of the first subband exciton. The re- 
sulting exciton-exciton resonance reduces the maximum ab- 
sorption strength and also shifts the optical transition energies. 
As an important consequence uncorrected photoluminescence 
overestimates the abundance of armchair-like tubes. We sug- 
gested several experiments to verify our ideas, among them 
infrared Raman spectroscopy in resonance with the first sub- 
band exciton. 
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FIG. 6: Photoluminescence map calculated without exciton-exciton resonance [B = in Eq.(l) of the paper]. The simulated sample is the 
same as described in the paper. This figure is to be compared with the calculation including exciton-exciton resonances [Fig. 4(a) of the main 
paper] and experiment [Bachilo et ah, Science 298, 2361 (2002)]. Note, in particular, the many peaks below the isotropic limit (dashed line), 
which are no longer visible when the eh22 — * 2eli22 decay is taken into account. 



TABLE I: Topological and optical parameters for single-walled carbon nanotubes. E\\ and E22 are the bare exciton energies as opposed to 
the measured renormalized transition energies hV22- The hV22 were taken from Ref. 2, 7, and 9 of the paper. PL int. (Raman int.) are the 
photoluminescence [ehn resonant Raman, see Ref. 21 of the paper] intensities; they are normalized to the intensity of the (5,4) nanotube. 
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